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Abstract. We compare the rate of kinetic energy injected by stellar winds into the extragalactic H ll region NGC 
2363 to the luminosity needed to feed the observed temperature fluctuations. The kinetic luminosity associated to 
the winds is estimated by means of two different evolutionary synthesis codes, one of which takes into account the 
statistical fluctuations expected in the Initial Mass Function. We find that, even in the most favorable conditions 
considered by our model, such luminosity is much smaller than the luminosity needed to account for the observed 
temperature fluctuations. The assumptions underlying our study are emphasized as possible sources of uncertainty 
affecting our results. 
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1. Introduction 

The presence of temperature fluctuations in photoionized 
regions has been a matter of debate since the pioneer- 
ing work by Peimbert (1967). Although little doubt on 
their relevance in real nebulae seems possible on observa- 



tional grounds (see the review by Peimbert 199f;), their 
inclusion in a theoretical framework is still controversial, 
given both the incompleteness of the present theoretical 
scenario, and the technical difficulties implied by their in- 
clusion in photoionization models. Plain photoionization 
theory predicts temperature fluctuations to be very small, 
mainly due to the steep dependence of the cooling rate 
on temperature, which implies that their lifetime would 
be quite short; yet, the existence of a mechanism steadily 
providing the energy required to feed them is not excluded 
by this line of reasoning. The importance of the question 
can be barely overrated given that, if temperature fluctu- 
ations turn out to be as big as observations suggest, the 
present determinations of chemical abundances should be 
rescaled by as much as -1-0.5 dex. 

For a long time, the only way temperature fluctuations 
were accounted for in theoretical models was acknowledg- 
ing the impossibility to reproduce the observed intensity of 



the most affected lines; see, e.g., Luridiana et al. (1999) for 
the case of NGC 2363, and Stasiiiska & Schaercr (199£) for 
the case of IZwl8. Recently, Binette & Luridiana (2000) 
developed a model to quantify the effect of temperature 
fluctuations on the ionization and energy balance of neb- 
ulae. In the present work, we use their schema to investi- 
gate whether the kinetic energy provided by stellar winds 
could feed the temperature fluctuations observed in NGC 
2363. Throughout the paper, by 'stellar winds' we generi- 
cally refer to all those phenomena involving the ejection of 
mass from stars into the interstellar medium (ISM). Given 
the age range considered, the only contributors to stellar 
winds actually included in the calculations are hot stars 
and Wolf-Rayet stars (WRs). 

The structure of the paper is as follows: in Sect. ||we 
summarize the relevant properties of NGC 2363. Sect. || 
describes the model used to represent temperature fluc- 
tuations in nebulae. In Sect. ^ we estimate the kinetic 
luminosity of the stellar cluster in NGC 2363, and apply 
the method of Sect. ^ to determine whether stellar winds 
could be responsible for the temperature fluctuations ob- 
served in NGC 2363. Finally, in Sect. ||we analyze several 
possible sources of uncertainty, and in Sect. ^ we summa- 
rize the main conclusions. 
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2. General properties of NGC 2363 

NGC 2363 is a giant cxtragalactic H ii region, located in 
the south-west end of the irregular galaxy NGC 2366. It 
is one of the brightest extragalactic H ii regions known, 
and an ideal subject for our study, since plenty of data 
are available, both observational and theoretical. Two dis- 
tinctive knots can be clearly distinguished in the region: 
throughout this paper, we will always refer to the bright- 
est, youngest knot, often called 'knot A'. 

NGC 2363 has been observed by several groups (e.g.. 



INUL^ nas been observed by several groups (e.g., 
Peimbert et al. 1986| : Gonzalez-Delgado et al. 1994|; Izotov 



et al. IQQ'il), and iu udL'k'd by [Lm'ldiaua eL al. (1999| ) and 



Drisseii cL al. (200C ) among uLhers, while Lhe LeniperaLure- 



fluctu ation parameter has been determined by jConzalez- 



Delgad o et al. (1994 ) in both knots 



In the present work, we will assume that both the stel- 
lar cluster and the gas nebula o f NGC 2363 are well de - 
scribed by the best-fit model of Luridiana et al. (1999| ), 
which is a spherical, hollow, radiation-bounded nebula, 
consisting of two concentric shells of different densities, 
ionized by a young cluster undergoing an extended burst 
of star formation. The properties of the model relevant for 
this work will be described throughout the text. 

3. The energy implications of temperature 
fluctuations 

To quantify the energy injection rate needed to fuel tem- 



peratu: 
propos 



fluctuationo in a nebula, wo follow the model 



d by B kietto fe Luridiana (2000 | ) and Binotto et 



al. (2001). The first of these two papers contains a de- 



tailed description of the method, while we refer to the 
second for an alternative representation of the results. The 
model depicts temperature fluctuations as a collection of 
hot spots, arising above a uniform equilibrium tempera- 
ture floor, and fed by an unknown heating agent. This 
representation allows one to compute the average thermal 
properties of the nebula in the temperature-fluctuation 
regime, but without having to compute real localized fluc- 
tuations, as described below. 

3.1. Definition of temperature fluctuations 

The amplitude of temperature fluctuations in a neb- 
ula is measured by the parameter t^, deflned as follows 
( [Peimbert 1967| ): 



JyiUv)-To)^N,ir)N,{r)dV 
T^J^N,ir)N,ir)dV 



(1) 



where Tc(r) is the local electron temperature, A'c(r) and 
A'i(r) are the local values of the electron and ionic density 
respectively, V is the observed volume, and the average 
temperature Tq is given by: 



Tn = 



/^re(r)jVe(r)jVi(r)dV^ 
JyN,ir)Ni{v)dV 



(2) 



Since the temperature of a model H ii region is not 
spatially constant, temperature fluctuations appear and 
evolve following the evolution of the ionization fleld even in 
the simplest case of a static, pure photoionization nebula; 



in particular, Perez (1997) showed that large temperature 
fluctuations arise during the WR phase, as a consequence 
of the hardening of the spectrum. Following the deflnition 
by Ferland (1996), we will refer to such 'structural' tem- 



perature fluctuations as t'^^^.. Typical t^^^. values for pho- 
toionization models of chemically and spatially homoge- 
neous model nebulae are in the range 0. 00 < ^^tr ^ 0.02. 



The b est-flt photoionization model by Luridiana et al 
(1999|) yields a value 4, - 0.009. 



On the observational side, Gonzalez-Delgado et al 
(1994D flnd for NGC 2363 the value t^^s = 0-064 by com- 



paring the Paschen temperature, Tc(Pa), to the [O ill] 
temperature, obtained from the A4363/A5007 ratio. Even 
if in this case the formal error on t'^^^ is quite large (we 
calculated cr(iobs) — 0-045, mainly due to the large er- 
ror on To (Pa)), the differences between the model predic- 
tions and the observations seem to imply that an extra- 
heating mechanism, other than photoionization, is at work 
in most objects, producing an additional i^xtra such that 
^2^^^ = ^2^^ + tl^^^^. In the case of NGC 2363, we find 



0.055 ±0.045. Note that such value still leaves the 



door open to the possibility, small but not negligible, that 
^oxtra — 0.00; however, we assume this not to be the case. 
In this respect, it is interesting to note that ^2^^ — 0.098 



in knot B of the same region (Gonzalez-Delgado et al 



1994 ) . Given the presumed age and metallicity of this re- 



gion ( Luridiana et al. 1999| ) , and assuming the same filling 
factor as in knot A, we can roughly estimate tl^^ = 0.02 
(see, e.g., [Perez 1997D , yielding tl^^^^ = tl^^-t^, = 0.078, 
with an associated error of about cr(iobs) ~ 0.019. 

More generally, even though the errors might be com- 
patible with ^oxtra = in individual cases of ionized re- 
gions with i2^j, > 0, this is certainly not the case when 
large samples are considered. As an example, galactic H ii 
regions and planetary nebulae collectively show higher 
values, typically around 0.04 (e.g., p^eimbert et al. 1995 ; 



Peimbert 1995, and references therein). 



The presence of temperature fiuctuations inside NGC 
2363 bears important cosmological consequences, as it af- 
fects the determination of the chemical abundances, low- 
ering the extrapolated primordial helium abundance Yp 
(e.g., Peimbert et al. 2001a, and references therein). In our 



case, <2j^g — 0.064 implies a downward change in the he- 
lium abundance of NGC 2363 of order 2%, the exact figure 
depending on which He lines the abundance determina- 
tion is based on (see also Gonzalez-Delgado et al. 1994). 
It is important to note, however, that temperature fiuctu- 
ations are not the only factor affecting a proper chemical 
abundance determination, other being, e.g, the ionization 
structure and the coUisional excitation of the hydrogen 



lines (Btasihska & Izotov 2001; Peimbert et al. 2001t) 
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3.2. The hot-spot scheme 

The overall energy balance in a nebula is described by the 
following equation: 



H = G 



(3) 



where H and G are the heating and cooling rates respec- 
tively. In a nebula in ionization equilibrium, with no extra- 
heating sources, the heating is provided by photoioniza- 
tion, and the cooling takes place mainly through collision- 
ally excited line emission and free-free radiation. A local 
equilibrium temperature Teq is implicitly defined by Eq. 
^, such that the heating terms counterbalance the cool- 
ing terms. This generic model is the reference tpxtra = 
model. 

In the hot-spot model, a particular temperature profile 
of a hypothetic nebula is postulated, consisting of a col- 
lection of randomly-generated hot spots arising above the 
Teq floor. Following the definitions of Eqs. |l| and ||, the 
^extra ^^'^ ^0 valucs Can be computed, giving i^xtra > 
and To > Teq. The recombination rates of this hypo- 
thetic model depend on the assumed temperature struc- 
ture, being in general different from those of the reference 
^cxtra = model. It is possible to account for such varia- 
tion by introducing a new 'global' temperature, Tj-oc, such 
that the intensity of a recombination line is given by: 




J-rcc '-^ J^roc 



(4) 



Fig. 1. Fhcat as a function of i^xtra ^'^^ the model of NGC 
2363 (solid line) and for a simpler model (dashed line, see 
text). The textra value inferred from observations of NGC 
2363 is also shown. 



^cq + Qcq is the equilibrium cooling rate, and Tfluc — -^eq 
gives the extra-luminosity radiated because of tempera- 
ture fluctuations, so that Fhcat gives a measure of the ex- 
cess energy provided by the putative heating mechanism 
driving the fluctuations, and dispersed through collisional- 
line enhancement. 

The relationship between i^xtra and Fhcat depends on 
the properties of both the ionizing source and the gaseous 
nebula. In Figur e [l| we show the depe ndence of Fhcat 
on igxtra fo'" the Luridiana et al. (1999| ) model of NGC 
2363, which we computed using a modified version of the 



^here g is a representative average of the H I recom- phptoionization-shock code MAPPINGS IC (|Ferruit et al 



binaticfrr 



range. 

Analogously, the intensity of a coUisionally excited line 
in the temperature-fluctuation regime can be calculated 
by means of a 'global' coUisional temperature Tcoi, such 
that the coUisional excitation and de-excitation rates are 



cxpoiieiiL <^(Tc) Hi Lhe approprlaLe LemperaLure ^^q^ ^ g^^^^^j ^^^^^ ^ ^^^^^^^ ^^^^1 ^^^^^^ 



proporlioual Lo T^l^ exp AE^j / kTcoi] and T^l^ I'espec- 
tively, Willi Lhe P and AEij values appropriate for each 
considered line. 



puted for each considered case, and the calibrations ob- 
tained for simpler models should only be used as rough 
guidelines when estimating the energy implications of tem- 
perature fluctuations in a given object. We can illustrate 
this point comparing our calibration to the one by Binette 



& Luridiana (2000) at a representative temperature. We 



Thiough Lhe defliiiliou of Lhese corrected Leiiipera- 
tures, Lhe new fine eiiilssivlLies caii be coiiipuLed iii a 
straightforward way. In the hot-spot model, the relation- 
ship between i^^tra and the derived temperatures is put 
in analytical form and generalized, so that the effects of 
fluctuations on the output quantities can be easily com- 
puted. For each postulated t^xtra value, a plain photoion- 
ization model is calculated, in which the line emissivities 



first introduce an equivalent effective temperature for the 
cluster, deflned, according to the method by Mas-Hesse 



Kunth (1991D, as the temperature of an early-type star 



are coir 

the global energy balance is consequently modified. This 
energetic change can be expressed by means of the quan- 
tity Fhcat, defined as follows: 



with the same Q(IIe")/(5(Il") ratio as the synthetic stel- 
lar energy dist ribution (SED). Using the calibration by 
Panagia (1973|) , we found an equivalent effective temper- 
ature Tcff = 43700 K. In Figure [l| we reproduce the cali- 
bration computed by Binette fc Luridiana (2000 ) for the 
case of a constant-density, Z=0.004 nebula io nized by an 
unblanketed L TE atmosphere of T=45,000 K (|Hummer fc 



ecLed for lIic Lemperamre-flucLuaLion ^necLS,^nd ^^-^^^^^^ Although the two curves show the same 



qualitative behavior, they rapidly diverge for increasing 



Tflu 



heat 



(5) 



where Tfluc and Lcq are the energies radiated by the 
nebula through line emission with and without temper- 
ature fluctuations respectively, and Qcq is the energy ra- 
diated through processes other than line emission. Thus, 



iextra ValueS. 



4. The kinetic luminosity of the stellar cluster 

In this section, we will compare the kinetic luminosity of 
NGC 2363 to the heating rate implied by the observed 
temperature fluctuations. To estimate the kinetic lumi- 
nosity produced by the stellar winds in NGC 2363, we 
considered a synthetic stellar population with the param- 
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Parameter 



Value 



Synthesis code IMF filling 



SFR 

1 + x 

Age (Myr) 
Z 

Mtot(lM0 



120 M 



Extended burst, At 

2.00 

120 

3.0 

0.004 

1.07 X 10^ Mq 



1.4 



Table 1. Adopted parameters for the stellar population 
in NGC 2363. 



eters resumed in Table [l]. The only differ ence with respect 
to the model by Luridiana et al. (1999| ) is the adoption 
oi Z = 0.004 instead oi Z — 0.005, to avoid interpola- 
tion between spectra, and to allow a more direct compar- 



ison with the models by Binette & Luridiana (2000). This 
choice has little influence on the temperature-fluctuation 
problem, since we estimate that adopting Z = 0.005 would 



increa se the kinetic luminosity by less than 50% ( Leithercr 



et al. ; 992), which is not sufficient to change, as we will 
see, the results of this work. 

Three different estimates of the kinetic luminosity L\^in 
and the total integrated wind energy for this source are 
listed in Table S. One of them was obtained with the first 



release of Starhurst99 (see Leithercr et al. 1999) run with a 
theoretical wind treatment, which is the default set by the 
code and the one preferred by the authors of Starhurst99. 
The other two were o btained with an updated ver sion of 
the synthesis cod e of Cervino fc Mas-Hesse (1994 ) (here- 
after CMHK, see ICervino et al. (2001a) for details), wit h 



the wind treatment illustrated in Cerviiio et al. (2001a ) 



The reasons for considering two different synthesis codes 
with two different mass-loss laws reside in the importance 
of showing the effects of different assumptions on the wind 
treatment (a possibility granted by the Starhurst99 code, 
which can be set to any of four different laws), together 
with the importance of estimating the expected statistical 
dispersion in the energy-related quantitities (a possibil- 
ity granted by the CMHK code, either analytically or by 
means of Monte Carlo simulations). Stated otherwise, it 
is not our intention to compare the two codes, but rather 
the two wind treatments: had we chosen the same one 



as Cerviiio et al. (2001a) in the Starhurst99 settings, we 
would have found the same result as with the CMHK code. 

The two estimates of Lkin obtained with the CMHK 
code have been obtained with an analytical and a Monte 
Carlo representation of the IMF respectively, and they 
both take into account the effect of statistical fluctuations 
in the IMF. When statistical effects in the IMF are taken 
into account, the population properties are no more univo- 
cally determined, being instead distributed along a prob- 
ability distribution curve. The characteristic parameters 
of the curve can be determined a nalytically or by mean s 
of Monte Carlo simulations (see Cerviiio et al. 2001b| ). 
The second line of Table || lists the analytical estimate 
of the average wind luminosity and the total integrated 
wind energy obtained with the CMHK code, as well as 
the corresponding 90% confldence intervals (i.e., the 5**^ 



(lO'^-' erg s- 



^'kin 
(10^2 



Starhurst99 

CMHK 

CMHK 



0.1 



Analytical 
Analytical 1 
Monte Carlo 1 



^^-0.23 
1f5+0.37 



5.0 

4.34^ 

A QO+0.8» 



0.68 
0.61 



Table 2. Kinetic luminosities and integrated kinetic en- 
ergies for the three considered SEDs. 



and the 95*"^ percentile). The third line quotes the same 
quantitites obtained by means of Monte Carlo simula- 
tions, again with the uncertainties corresponding to the 
90% confidence interval. These two estimates of Lkin are 
in excellent agreement, showing the consistency of the two 
statistical approaches of the CMHK code. In all the cases, 
the extended burst scenario has been represented as the 
sum of individual instantaneous bursts of different ages 
and roughly equal masses. In the case of the Monte Carlo 
method, different simulations, taken from an ensemble of 
5000 runs of 10^ stars each, have been summed up until 
the mass of each individual cluster was reached. The to- 
tal number of independent Monte Carlo simulations of the 
extended burst obtained through this procedure was 64. 

The average estimates obtained with CMHK give 
roughly 70% higher luminosity than the Starburst99 value; 
this is mainly due to the difference between the wind treat- 
ment adopted in the two codes, and is not directly related 
to differences in the spectra, which are, for the scopes of 
the present work, largely negligible. The integrated en- 
ergy value obtained with Starburst99 is higher than the 
CMHK one, a somewhat surprising trend given the corre- 
sponding luminosity values. These figures depend on the 
winds evolution with time: the wind treatment adopted in 
the Starburst99 calculation gives less energetic winds than 
the one implemented in the CMHK code at earlier ages, 
while for ages greater than about 2.5 Myr the relationship 
is inverted. 

It is now possible to compute the Fheat values associ- 
ated with the kinetic energy input, according to the gen- 
eral expression: 



pkin 
heat 



Lk 



(6) 



For simplicity, only one value for the equilibrium cool- 
ing rate will be used in t he T^"]^t comp utations, resulting 
from a CLOUDY model ( Ferland 1996| ) of the region run 
with the Starburst99 SED: 



Leq + Qcq - G~s = 1.46 X lO^ierg s" 



(7) 



and only the amount of kinetic luminosity deposited ac- 
cording to each SED will be varied. This introduces only 
a minor approximation, since the differences between syn- 
thesis codes, as well as the stochasticity of the IMF, affect 
the shape of the spectrum only marginally before the WR 
phase, which, in our extended-burst model accounts only 
for a minor fraction of the flux. This statement is not valid 
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Synthesis code IMF fliling 



- heat 



tt„ X 10^ 



Starburst99 

CMHK 

CMHK 



Analytical 
Analytical 
Monte Carlo 



0.0047 0.60 

omietoMil 0.96: 

0.0076^o^;;i? 0.96: 



0.19 
0.40 



Table 3. T^^"l^ values associated to the wind luminosities, 



and inferred t^j^^ values. 



An independent estimate of the efficiency of the ther- 
malization of the wind energy can be done by considering 
the X-ray component in the spectrum of NGC 2363 de- 
tected by ^tevens fc Strickland (1998| ) with ROSAT; they 
found Lx > 6.6 x 10'^'^ erg s~^ assuming a distance of 3.44 
Mpc to the object, which becomes ^ 7.3 x 10^^ erg s 



rescaling to the distance of 3.8 Mpc assumed by Luridiana 



et al. (1999). Assuming that the origin of the X-ray com- 



in general, since for WR-dominated bursts, a large disper- 
sion is e xpected in the ionizing flux, hence in the cooling 
rate: 



see 
The 



Cervifio et al. (200 IbD 



ponent is the reprocessing of the wind kinetic energy, we 
can infer typical values for the thermalization efficiency of 
the order > 10%. 



^heat values obtained are listed in Table ^, to- 



gether with the inferred t^^^^^ obtained from the hot-spot 5.2. Pure photoionization component 



J); we will refer to these 
to emphasize that they 



model of NGC 2363 (cf. Figure 
^extra values with the symbol t^j, 
are inferred under the assumption that the temperature 
fluctuations are driven by the stellar-wind kinetic energy. 
All the t^jjj values obtained are extremely small, well below 
the tgxtra inferred from observations. However, there are 
still a few issues to consider before drawing any conclusion 
from the comparison between the t^^^ values of Table |^ and 
the t^^tra value of NGC 2363 inferred in Sect. O. 



5. Discussion 

5.1. Wind-luminosity tfiermalization 

An important issue to consider is the thermalization ef- 
ficiency of the wind kinetic energy. Even assuming that 
our estimates of the wind luminosities in NGC 2363 are 
highly accurate, the derived i^jj^ values should be corrected 
downward to take into account the fraction of the energy 
injected to the ISM which is eventually thermalized. 

1-D hydrodynamical models of the bubble around OB 
associations indicate that such fraction is about 80% 



(Pliischke 2001; Phischke et al., in preparation). An es- 



timate of the thermalization efficiency in the case of NGC 
2363 can be made by comparing the total kinetic energy 
injected into the region since the begi nning of the star - 
burst to the observed kinetic energy. Roy et al. (1991 ) 
observed in NGC 2363 a bubble with a 200 pc diameter, 
expanding with a velocity of 45 km s~^, and calculated 
that the total kinetic energy involved is E°}' ^ ~ 2 x 10 ^^ 



ergs, in agreement with the estimate made by Luridiana et 



al. (1999 ) . Roy et al. (1992 ) found a high- velocity compo- 
nent in the gas, with a kinetic energy of E^f^ ^ 10^^ ergs. 



Gonzalez-Delgado et al. (1994 ) confirmed the presence of 



this high-velocity gas, and estimated a kinetic energy of 
^khT ~ 3 X 10^2 grgg^ 

These figures should be compared to the estimated val- 
ues of the integrated wind kinetic energy, which are listed 
in Table || for the three considered cases. If we add the 
kinetic ener gy of the low- velocity bubble observed by Roy 



et al. (1991|) to that of the hi gh-velocity gas observed by 



Gonzalez-Delgado ct al. (1994| ), we find that the efficiency 



of the thermalization of wind kinetic energy is close to 0, 
lowering further the computed f^j^^ values. 



The value of i^^j. could be higher if the ionizing spectrum 
turned out to be harder than supposed, as implied by 
the X-ray component detected by ^tevens fc Stricklanc^ 
|(1998| ) ; in this case, the need for an extra-heating source to 
be added to the photoionization model would be propor- 
tionally smaller. To investigate this possibility, we com- 
puted a Cloudy model with the analytical SED of the 
CMHK code, modified to account for the transformation 
of 20% of the kinetic energy into X rays ( |Cervino et al 



in preparation). This experimental model is energetically 
equivalent to a hot-spot model in which the extra heat- 
ing is provided by the wind kinetic luminosity thermalized 
with a 20% efficiency. Indeed, we found that for this model 



ig^j. = 0.009, i.e. the same as the Luridiana et al. (1999) 
to within 0.001, confirming that the wind kinetic luminos- 
ity is energetically insufficient to account for the observed 
temperature fluctuations. 

5.3. Influence of stellar rotation 



As it has been pointed out by Meynet & Maeder (2000), 
rotation dramatically changes the properties of massive- 
star models. In particular, it increases the mechanical en- 
ergy released to the ISM (see, e.g., Maeder fc Meynet 
^001 ). Additionally, the winds of rotating massive stars 
can be highly non-spherical, with strong polar or equa- 
torial structures, depending on the effective temperature 
and angular velocity of the star (Maeder & Desjacques 

pooi| ). 

Rotation makes it possible to explain in several cases 
the observed features of individual stars, but, unfortu- 
nately, the evaluation of its effects on the integrated spec- 
trum of star forming regions is not feasible yet. First, the 
distribution of angular velocities of the stars in the cluster 
should be known. Second, the distribution of inclination 
angles should also be known, since the emitted luminosity 
depends on the inclination angle. Third, it would be neces- 
sary to establish new homology relations in order to obtain 
the isochrones needed to compute the emission spectrum 
at any given time. 

Thus, we can only predict qualitatively that rotating 
stars certainly produce a turbulent interstellar medium, 
possibly with strong density and temperature inhomo- 
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geneities, due to both the increased wind energy and the 
anisotropy of their winds. From the point of view of our 
study, the increased energy injected into the medium, 
would translate into a correspondingly higher r}jj,"j value, 
and the turbulence created in the ISM would possibly in- 
crease the thermalization efficiency of such energy, in such 
a way that the 'actual' Fhcat of the region would approach 



its upper limit Fj^'J^^^. (see also Sect. 5.1). 

Thus, at a qualitative level, the effect of stellar tracks 
with rotation on our temperature- fluctuation model would 
be to increase the temperature fluctuation amplitude the- 
oretically achievable through energy injection by stellar 
winds, as compared to non-rotating models. 

5.4. Distance to NGC 2366 

The distance to NGC 2366 plays several roles in our anal- 
ysis. First, the average properties of a photoionization 
model constrained by observational data depend on the 
assumed distance in complex ways. Second, a smaller dis- 
tance implies a smaller region, hence a relatively larger 
statistical dispersion. Third, the hot-spot model is cali- 
brated to a specific model, so that, should the photoioniza- 
tion model change, the hot-spot model would also change. 

In our analysis, based on the photoionization model 
by Luridiana et al. (1999), we assumed for the dis- 



tance the value 3.8 Mpc determined by Sandage fc 



Tammann (197(: ) . However, there are a number of more re- 
cent distance determinations indicating smaller distances: 
Tikhonov et al. (1991 ) derived a distance of 3.4 Mpc to 
NGC 236 6 through photograph ic photometry of its bright- 
est stars; Aparicio et al. (1995 ) obtained the valu e of 2.9 



Mpc with CC D photometry of its brightest stars; Tolstoy 
et al. (1995 ) determined a dis tance of 3.44 M pc from 

obtained 



Cepheid light curves and colors; Jurcevic 
the value 3.73 Mpc, with an associated error of 0.04 dex, 
from a study of the period-luminosity relationship of the 
red supergiant variables in NGC 2366. 

To illustrate how our conclusions would change un- 
der a different assumption on the distance, we consider a 
value of 3.4 Mpc, i.e. ^--^10% less than the assumed dis- 
tance. Given the observed II/3 intensity, the smaller dis- 
tance yields a 20% smaller rate of ionizing photons, Q(II'^). 
Since in this case the region turns out to be smaller both 
in linear and in angular dimensions, more input param- 
eters should be changed, in order to fulfill the observa- 
tional constraint on the observed size of the nebula, even 
if the remaining constraints, such as the relevant line ra- 
tios, change negligibly. Though it is beyond the scope 
of this paper to calculate a full revised photoionization 
model, we can expect that a satisfactory fit could be ob- 
tained by means of relatively small adjustments in the 
density structure of the nebula surrounding the ionizing 
cluster, with the model stellar population rescaled to a to- 
tal mass MjQj = 0.80 x Mtot, and the other stellar param- 
eters (e.g., SFR, IMF, etc.) left unchanged. These changes 
would translate into a F}^™^ smaller by 20%, with a sta- 



tistical dispersion proportionally larger, due to the larger 
relative weight of the statistical fluctuations in the cluster. 

The hot-spot model should also be accordingly mod- 
ified, to take into account the properties of the revised 
model. However, we don't expect it too change too much, 
since the ionization parameter of the revised model would 
be essentially the same as the one of the old model. 
Summarizing, we estimate that a downward revision of 
the assumed distance would not significantly change our 
conclusions. 

5.5. Temperature-fluctuation profile 

In the interpretation of our results, it is important to take 
into account that they were obtained for a rather specific 
temperature-fluctuation pattern. As stated by Binctte & 



Luridiana (2000), the model returns the correct results for 



fluctuations resembling those depicted in their Figure 1, 
but for radically different patterns of hot spots (different 
in frequency, width, and/or amplitudes) the model would 
only provide a first order estimate of the relationship be- 
tween and Fhcat; Binette & Luridiana (200C) estimate 



that the total uncertainty on Fhcat resulting from this ap- 
proximation is less than 20%. 

6. Conclusions 

Our results suggest that the kinetic energy provided to the 
ISM by the stars through stellar winds cannot account for 
the observed temperature fluctuations in NGC 2363. This 
result holds even if a thermalization efficiency of 100% is 
assumed; however, the comparison of the observed kinetic 
energy with the theoretical estimate of the wind kinetic 
energy suggest that such efficiency is rather low. These 



results confirm the conclusion drawn by Binctte et al 



(2001 ), and leave the question of the nature of energy 



source of temperature fluctuations open. New insights into 
the problem could possibly come from the use of stel- 
lar tracks with rotation, and/or the consideration of a 
temperature-fluctuation pattern radically different from 
the one used in the present work. 
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